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Abstract 

We analyze various limits of vibrationally coupled resonant electron transport in single-molecule 
junctions. Based on a master equation approach, we discuss analytic and numerical results for 
junctions under a high bias voltage or weak electronic-vibrational coupling. It is shown that in 
these limits the vibrational excitation of the molecular bridge increases indefinitely, i.e. the junction 
exhibits a vibrational instability. Moreover, our analysis provides analytic results for the vibrational 
distribution function and reveals that these vibrational instabilities are related to electron-hole pair 
creation processes. 

PACS numbers: 73.23.-b,85.65.+h,71.38.-k 
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I. INTRODUCTION 



Charge transport through nanostructures has been of great interest ever since nanofabri- 
cation techniques have been emerged jl-9]. The quantum mechanical nature of the charge- 
carriers in these nanostructures gives rise to many intriguing transport phenomena, e.g. 
strongly nonlinear transport characteristics {4, 5|. When these structures became contin- 
uously smaller, it was realized that vibrational (or phononic) de gree s of freedom play an 



ft 

of nanoelectronics is found in molecular electronic devices |7H9, 



important role in this nonequilibrium transport problem 



Of. The ultimate limit 
12J, where the building 



blocks consist of single molecules. Due to their small size and mass, molecules often show 



strong correlations between their electronic and vibrational degrees of freedom 



13|-|25|. It 



is thus of great interest to understand electron transport through a nanostructure that ex- 



hibits elec 
electrode 



;ronic-vibrational coupling, such as a single molecule coupled to a left and a right 



12 



26j . Analyzing the limits of this transport problem, especially in the resonant 



transport regime, facilitates the understanding of nonequilibrium transport at the nanoscale. 

Due to electronic- vibrational coupling, transport through a single-molecule junction com- 
prises not only charge-exchange but also energy-exchange processes with the leads. Examples 
of such processes are depicted in Fig. [TJ Thereby, Panel a) and b) represent transport pro- 
cesses, where the tunneling electron vibrationally excites (heating) and deexcites (cooling) 
the molecular bridge, respectively. Panel c) shows a process, where in two sequential tunnel- 
ing processes an electron-hole pair is created in the left lead upon absorption of vibrational 
energy from the molecular bridge. Such energy-exchange processes with the leads play an 
important role in molecular junctions. This can be illustrated by considering a junction, 
where the temperature in the leads, T, represents the largest energy scale. In this limit the 
population of the molecular energy levels is determined by the thermal distribution that also 
determines the population of the levels in the leads. The corresponding average vibrational 
excitation of the molecular bridge thus increases indefinitely as T — > 00. That way, the 
high-temperature limit exhibits a vibrational instability in a trivial sense. Similarly, the 
static limit [27] , where the frequency of the vibrational modes represent the smallest energy 
scale, Q — > {k^T 3> Q), results in an infinite vibrational excitation, since energy-exchange 
processes transfer the thermal excitation of the leads to the vibrational degrees of freedom 
of the molecular bridge. 
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FIG. 1: Basic processes in vibrationally coupled electron transport through single-molecule junc- 
tions. Panel a) (Panel b)) depicts an example for an emission (absorption) process, where an 
electron sequentially tunnels from the left lead onto the molecule and further to the right lead, 
thereby singly exciting (deexciting) the vibrational mode of the molecular bridge (red (blue) wiggly 
line). The frequency of the vibrational mode is denoted by f2. An example for a process, where 
an electron-hole pair is created in the left lead by singly deexciting the vibrational mode of the 
bridge, is shown in Panel c). 



In this article, we focus on the limit of a large bias voltage, $ — > oo, as well as a weak 
electronic-vibrational coupling, A — > 0, and analyze the vibrational instabilities inherent to 
these limits. For an infinite bias voltage we find the corresponding vibrational excitation 
to diverge, because electron-hole pair creation processes are completely suppressed. This 
important cooling mechanism is also missing for vanishing electronic-vibrational coupling, 
A —7- 0, if the bias voltage is large enough. As a result, the limit of weak electronic-vibrational 
coupling exhibits another vibrational instability. This intriguing phenomenon was already 
reported and analyzed by Koch et al. 
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. Resonant absorption processes with a higher- 



lying electronic state may prevent this instability [29[. For moderate bias voltages that allow 
for the leading-order electron- hole pair creation processes in the limit A — > 0, we derive an 
analytic expression for the respective vibrational distribution function, which shows that the 
corresponding average vibrational excitation is finite and nonzero. Only in the off-resonant 
transport regime, i.e. for small bias voltages, we find that vibrational excitation vanishes as 
A -> 0. 



To study these limits we employ a master equation formalism 30H42j that is based on a 



second-order expansion in the coupling of the molecular bridge to the leads. This allows a 
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description of all resonant transport processes, but misses higher-order processes such as the 
broadening of levels due to the coupling between the molecule and the leads or equivalently 
co-tunneling processes 



43 



46| . Such processes can, in principle, be described by approaches 



that take into account higher-order effects, including advanced master equation approaches 



30, 



33 



methods 



4J42L; 
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43, scattering theory approaches 48M51|. nonequilibrium Green's function 



45 
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59j . or numerically exact methodologies 60H64j . However, since we 



are primarily interested in the limits of resonant electron transport through single-molecule 
junctions, higher-order processes play no role for our considerations. 

The article is organized as follows. In Section [III we introduce the Hamiltonian (Sec. Ill Ap 
and the master equation approach (Sec. Ill Bh that we use to describe vibrationally coupled 
electron transport through single-molecule junctions in the resonant transport regime. We 
present our results in Sec. IHIt wherein the transport characteristics of a molecular junction 
are analyzed in the limit of an infinite bias voltage (Sec. IIII Aj) and for a vanishing electronic- 
vibrational coupling (Sec. IIIIBl) . The latter case is studied in three different regimes: in Sec. 
IIII B II for a high bias voltage, where electron-hole pair creation processes are suppressed, in 
Sec. IIII B 21 for a lower bias voltage that allows for electron-hole pair creation processes, and 
in Sec. IIII B 31 including a second higher-lying electronic state that is vibrationally coupled. 



II. THEORY 



A. Model Hamiltonian 



We investigate vibrationally coupled electron transport through single- molecule junctions 
using the following Hamiltonian (throughout the article we use units where h = 1) 

H = ^e;cjc; + ^ e kAc k + ^(V fei cJ.Cj + ll.C.) 
i&A fceL,R k,i 

+fia t a + ^A i (a + a t ) (c\a - 8A , (1) 

where electronic states with energies 6j, located at the molecular bridge (M), are coupled 
by interaction matrix elements to electronic states in the leads (L,R). The energies 
of the lead states are labelled by e^. The operators cj/cj (cj,/cfc) denote respective cre- 
ation/annihilation operators for the states of the molecular bridge (leads). To model the 
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vibrational degrees of freedom of the junction we consider a single harmonic mode with fre- 
quency Q that is described by creation/annihilation operators a' /a. The coupling strength 
between the harmonic mode and the ith state of the molecular bridge is denoted by Aj. 
Thereby, the fixed parameters 5i ensure that there is no electronic-vibrational coupling in 



the electronic ground-state of the neutral molecular junction 26|, |65|, |66). This means that 



5i = (Si = 1), if the ith electronic state is unoccupied (occupied) in this reference state. 



We prediagonalize H by the small polaron transformation 
resulting in 

A; A 



Z$,H^H = Hs+H B + H : 



SB, 



i i<j 

Hb = ^e k c\c hl (3) 

k 

H S b = ^(VkX4c, + h.a). (4) 

ki 

Thereby, H s comprises the degrees of freedom of the molecular bridge (S), in particular the 
polaron-shifted electronic states, ^ = q + (25j — l)(Af/f2), vibrationally induced electron- 
electron interactions, ~ 2AjAj /0, and the harmonic mode. The leads degrees of freedom (B) 
are summarized in Hb- The term Hsb describes the coupling between the molecule and the 
leads, which, as a result of electronic-vibrational coupling, is renormalized by shift operators 
Xi = exp((Xi/Q)(a — a^)). Due to the small polaron transformation, there is no explicit 
electronic-vibrational coupling term in H$. 



B. Transport Theory 



The current and vibrational excitation of a single- molecule junction are obtained from 
the reduced density matrix p, which is given as the stationary limit of the well-established 



equation of motion 

dp® 
dt 

with 



26 



29 



32 



34 



67 



68| 



POO 

[H s ,p(t)] - / drtr B {[F SB , [H SB (T),p(t)p B ]]}, 
Jo 



H sb (t) = e - i{Hs+HB)T H S Be i{Hs+HB)T . 



(5) 



(6) 



Here, pe represents the equilibrium density matrix of the leads. Eq. ([5]) can be derived 



from the Nakajima-Zwanzig equation 
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70j, employing a second-order expansion in the 



coupling Hsb along with the so-called Markov approximation 

p(t - r) « e msT p(t)e- t71sT . (7) 

The master equation, Eq. (EJ, is evaluated assuming a stationary state, ^ = 0, and with 
basis functions \a)\v) that span the subspace of the electronic \a) and the vibrational degrees 
of freedom respectively. Thereby, the electronic basis functions are given in the occupa- 
tion number representation i.e. \a) = |ra 1 n 2 ..), where rij G {0, 1} denotes the population of 
the ith electronic state. The vibrational basis function \u) with v e No represents the z/th 
level of the harmonic mode. The coefficients of the reduced density matrix are thus denoted 
by 

PZ? = (*\P V1V2 W) = (a\(vMv2)\a'). (8) 

In the evaluation of Eq. (jSJ), we neglect principal value terms that describe the renormal- 
ization of the molecular energy levels due to the coupling between the bridge and the leads 



2y, |34|. These contributions are irrelevant for the results discussed in this work. 

Having determined the coefficients of the reduced density matrix, we can readily obtain 
the vibrational distribution function 



Pv 



The corresponding average vibrational excitation for a molecular junction with a single 
electronic state is given by 

A 2 

(a ] a) H = {a f a)jj + —(n 1 -2S 1 n 1 + S 1 ), (10) 
A 2 

= yVpi; + ™ ("l - 25ini + Si), 

V 

with 



(c\ci)h = J2pZ- (11) 



For transport through a molecular junction with two electronic states, the average vibra- 
tional excitation is calculated according to 

A 2 A 2 
i)h = Y^ vPv + ~ 25lHl + 5 ^ + ftf (™ 2 ~ 252712 + ^ ( 12 ) 

IF 

n 2 



+2^^(^/011,11 - hni - 6 x n 2 + 5i5 2 ) 



with 

ni = (c\c x ) H = ^2pu,u +Pw,w, (13) 
n 2 = (4c 2 ) H = J^Pn.n +Poi,oi- ( 14 ) 
The current through the molecular junction is obtained from the formula 

/■oo 

I = -i drtr s+B { [#sb (r), pp B ] /}, (15) 
Jo 

with / = — 2e^ SfcgL c l c fc- For derivation and evaluation of Eq. ( !T5|) . the same approx- 
imations as for Eq. (j5| have been used. 

III. RESULTS 

In this section, we use the methodology outlined in Sec. Ill Bl to investigate various limits 
of resonant electron transport through single molecules. In particular, we study the current 
and the vibrational excitation of a molecular junction in the limit of an infinite bias voltage 
(Sec. Ill I A I) and for vanishing vibronic coupling (Sec. lIHB) . Our analysis includes analytic as 



well as numerical results. For the latter we have used a set of 400 vibrational basis functions 
to obtain numerical convergence. Since we consider systems without quasi-degeneracies, 
coherences of the density matrix are negligible, i.e. p v ^ v J = p^^^i^^a' ■ Furthermore, we 
use the wide-band approximation and assume the bias voltage $ to drop symmetrically at 
the contacts. 

A. The high-bias limit — > oo 

In the resonant transport regime, the average vibrational excitation of a single-molecule 



junction typically increases very rapidly with increasing bias voltage [26|, |29|, |45| . An example 
for this behavior is given in Fig. |2k., which shows the vibrational excitation of a molecular 
junction induced by inelastic transport processes through a single electronic state (like the 
one depicted in Fig. Hk). Increasing the bias voltage, more and more inelastic transport 
processes become active. However, as they also involve an increasing number of vibrational 
quanta, which results in an unfavorable Franck-Condon overlap, they are typically strongly 
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suppressed. Hence, the question arises wether the level of vibrational excitation, induced by 
these processes, does or does not saturate in the limit $ — > oo. 

To prove that vibrational excitation increases indefinitely in the limit $ — > oo, we as- 
sume in the following that vibrational excitation is finite and derive a contradiction to this 
assumption. To this end, we evaluate Eq. (j5J) between the basis functions (0|(^i| and |i/i)|0), 

EM^ + n(u 3 - vi))t k \x 1>v1V3 \ 2 p ™ = ( 16 ) 

5> - f K (ex + n(i* - ViWk^^p^. 

Multiplying this equation by v\, and summing up all equations, we obtain the following 
equation 

E vif K {ex + £l(v 3 - v^TxlX^fp^ = (17) 

Kv\v A 

If we assume vibrational excitation (a^a) to converge in the limit $ — > oo, so must the Ihs 
and the rhs of the above equation. Hence, we can take the limit $ — > oo on both sides in 
Eq. ifPTl) . replacing by 1 and /r(c) by 0, which results in 

^r L \x 1)U1U3 1 2 Polo 1 = v 1 r R \x 1:Ulua \ 2 p l [ 3 j 3 . (is) 

Applying the sum rule \X U ^\ 2 = 1 to the Ihs, and X^H^/J 2 = A 2 /S7 2 + /i to the r/is 
gives 

E ^ r ^oT = E^ + a 2 /^)^ 3 . (19) 

Analogously, using the (l|(z/i|...|i/i)|l)-projection of Eq. (|5]), we obtain: 

E(^ + A 2 /fi 2 )r L PoT = E ^ r RP"r- (2°) 

Subtracting Eq. (120]) from Eq. f fl9|) leads to the following equation 

- \ 2 /vt 2 E r L PoT = A V^ 2 E r ^T. ( 21 ) 

^3 ^3 

which is a contradiction, since the //is of Eq. (12T1) is negative while its rhs is positive, 
if A 7^ 0. Hence, for finite electronic- vibrational coupling A, vibrational excitation must 
diverge as $ — > oo. 
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Fig. [2] shows numerical results that illustrate this analytic finding. For these calculations, 
we have employed a model for a molecular junction that consists of a single electronic state 
at si = 0.15 eV and a single vibrational mode with frequency Q = 0.1 eV. The corresponding 
electronic-vibrational coupling strength is A = 0.06 eV. These parameters represent typical 
values for molecular junctions similar to those that have been employed in first-principles 
T^. The level-width functions, T L ,fL,r R = 67/zeV, are chosen to be much 



models 



66 



71 



smaller than the thermal broadening k^T = 1 meV [77|]. The black, gray and red bars in Fig. 
[2b represent the corresponding vibrational distribution function for bias voltages $ = 0.75 V, 
1.25 V and 1.75 V, respectively. It is seen that the vibrational distribution function becomes 
broader with increasing bias voltage and exhibits equal population of the lower vibrational 
states for large bias voltages. This leads to a continuous increase in vibrational excitation 
with increasing bias voltage, which does not saturate. 

We attribute this behavior to the lack of electron-hole pair creation processes (cf. Fig. [lb), 
which for Q ^> k^T can only deexcite the vibrational mode 26|, and which are completely 
blocked in the high bias limit. As a consequence, the number of excitation and deexcitation 
processes upon electron transport through the molecule is equal. The respective stationary 
state is given by a vibrational distribution function, where all vibrational levels are equally 
populated (cf. Fig. |2b). This can be rationalized considering the population of the z/th vibra- 
tional level p u after any transport process in this limit. Since there are as many excitation as 
deexcitation processes, the respective population is given by the sum Y^>=iVv'\Xvv'\ 2 with 
ft u the corresponding population before such a transport process. The only nonequilibrium 
state, which is invariant under these conditions, i.e. p u = p u , is that, where all vibrational 
levels are equally occupied since Y^=i \X W '\ 2 = 1- 

The respective current is given by 2erLrR/(rL + rR), which is the same result as Gurvitz 



et al. 



found for a junction in the high-bias limit without electronic-vibrational coupling. 



B. Weak electronic-vibrational coupling A — > 



For a finite bias voltage that fulfills e$ > 2(ei + Q), the vibrational excitation may 
also diverge in the limit of vanishing yibronic coupling A — > 0. This counter- intuitive phe- 



nomenon was reported before 



28 



3l| . In subsection IHIB 1[ we give a short overview of 



the phenomenon, and reinterpret the phenomenon in terms of electron-hole pair creation 
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FIG. 2: (Color online) Upper Panel: Average vibrational excitation, (a^a), of a molecular junction 
with a single electronic state that is coupled to a single vibrational mode. The inset shows the 
respective current-voltage characteristics. Lower Panel: Population of vibrational levels, p u , corre- 
sponding to three different values of the bias voltage Increasing the bias voltage, the vibrational 
distribution function becomes broader, and thus, the level of vibrational excitation gets higher. In 
the limit — > oo this leads to an infinite vibrational excitation. 

processes. This way, we establish the relation between the limit of a weak vibronic coupling 
and the high-bias limit discussed in the previous section. In the two subsequent sections, 
we investigate the limit A — > for lower bias voltages (Sec. IIII B 2p . where the leading- 
order electron-hole pair creation processes are not blocked, and in the presence of a second 
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higher-lying electronic state that is vibrationally coupled (Sec. IIHB 3j) . 



1. Without cooling by electron-hole pair creation processes 

The solid black line in Fig. [3] illustrates the phenomenon of vibrational instability in the 
limit A 0. It represents the vibrational excitation, resulting from vibrationally coupled 
transport through a single electronic state, as a function of the electronic- vibrational coupling 
strength A at a fixed bias voltage e$ = 0.55 eV > 2(ei + Vt). The model parameters are the 

Upon reducing A, the level of vibrational excitation 
where the slope of this increase gets larger for smaller 



same as those described in Sec. II 



monotonously increases [26 



28 



29 



A 



values of A. Moreover, the vibrational distribution function (data not shown) becomes 
broader and shows an increasing number of vibrational levels, starting from the ground-state, 
that are equally populated. Interestingly, this phenomenon occurs only if the associated bias 
voltage e$ exceeds 2(?i + Q). For such bias voltages electron-hole pair creation processes 
(cf. Fig. [lb) are suppressed to lowest order in A. In this sense, the limit A — > is equivalent 
to the high-bias limit, where also all relevant electron- hole pair creation processes become 
blocked and the level of vibrational excitation diverges for $ — > oo. 

These numerical findings can only give an example of the phenomenon, and are, further- 
more, limited by the number of basis functions employed. In the following, we therefore 
analyze this behavior for A — > in an analytic and more general way. To second order in A, 
with e$ > 2(ei + Q) > 2Q, we obtain from Eq. (|5]) the following set of equations [y > 1) 

o = r LP - - r RP - - r R {yX 2 P \f»- x - (2* + i)\ 2 P z + iy + i)a 2 p5T +1 ) > ( 22 ) 
o = r LP - - r RP - + r L (ux'dff"- 1 - (2» + i)x 2 P ^ + {v + i)\ 2 P ^ +1 ) . (23) 

From the latter equations one infers that the difference TlPq q — ^k P i\ ls °f second order 
in A. Therefore, we can replace the terms Y-^ P W by T^ P q v q in Eq. f l22|) (or vice versa in Eq. 
( 123]) ). Subtracting Eq. ( 122|) from Eq. ( 123|) (or vice versa) thus gives 

= (21/ + l) P Z ~ {v + 1KT +1 - (24) 

where a G {0, 1}. The recurrence relation defined by Eq. ( 1241) leads to divergent populations 
P v a u a ±00, if p^ a — P ®° a 7^ 0. Thus, the only solution, which is normalizable, is the one, 
where all vibrational levels are equally occupied, P a\ lu+1 ~ f> V a a = 0- This corresponds to an 
infinite vibrational excitation or a vibrational instability in the limit A — > 0. 
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FIG. 3: (Color online) Vibrational excitation as a function of the electronic-vibrational coupling 
strength A. The gray line shows the vibrational excitation obtained for the model system with a 
single electronic state at a bias voltage 2(ei + fi) > e$ > 2ei, reflecting the result of Sec. HUB 21 
The black line represents the phenomenon of vibrational instability, which is outlined in Sec. lIIIrTTl 
The dashed blue line shows that a higher-lying electronic state, which is coupled to the vibrational 
mode, prevents the vibrational instability depicted by the black line (cf. Sec. HUB 3|) . 

It is noted that rigorous divergence of the vibrational excitation is only found for an 
isolated molecular vibration, as it is described by the standard model of vibrationally cou- 
pled electron transport in molecular junctions considered here. In real molecular junctions, 
vibrational relaxation processes, introduced e.g. by coupling to phonons of the electrodes 
or other vibrational modes, would restrict the vibrational excitation to a finite value. How- 
ever, even in the presence of such relaxation mechanisms, vibrational excitation may not 
only monotonously increase with A, but may also decrease with an increasing electronic- 
vibrational coupling strength. 

As in the high-bias limit, the current obtained for A — > is also given by 2eFLrR/(rL + 
Tr). For weak electronic-vibrational coupling, vibrational processes (Fig. [I]) do not con- 
tribute to the current, as they take place on time scales much longer than electronic transport 
processes, i.e. processes that do not include an energy-exchange of the traversing electron 
with the vibrational mode of the molecular bridge. 
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2. With cooling by electron-hole pair creation processes 



For bias voltages in the range 2(ei + Q) > e$ > 2e±, the vibrational excitation of a mode 
coupled to a single electronic state remains finite in the limit A — > 0. This is illustrated 
by the gray line in Fig. El which has been obtained for the same parameters as considered 
above (black line) but with a smaller bias voltage $ = 0.45 V. Intriguingly, for A — > 0, the 
level of vibrational excitation approaches a non-zero value. 

This result can be derived from the master equation dSJ). For voltages with 2(ei + Q) > 
e$ > 2ei > Q, we obtain the following set of equations 

= r L (z/ 1 + l)A 2 poT + fa + l)A 2 r L ^+^ +1 (25) 

+r R (z/ 1 + i)x 2 p^ +1 - r R (2 I / 1 + i)x 2 P ^ + rViAVft- 1 " 1 - 1 , 

a- = r L ^Avr - r L (-(2^1 + i)A 2 Por + fa + i)a 2 po: +1 ^ i+1 ) . 

with 

A^ = rLP^-rR^r. (26) 

From Eqs. ( |25l) . which are valid to second order in A, we further deduce 

o = (u, (2 + r L /r R ) + v x + 1) p£* (27) 
-fa + i)A 2 (r L /r R + 2) - ^A^" 1 , 

where a6 {0,1}. The solution of this equation is given by 

pZ = p^V (2 + r L /r R ) , (28) 

pu = pi:, + p\i (29) 

r L + r R ( r R 



V 

) 



r L + 2r R vr L + 2r R/ 

which corresponds to an average vibrational excitation of 

(a+a) = r R / (r L + r R ) r ^ R 1/2. (30) 

Thus, vibrational excitation approaches a finite value as A — > and does not vanish, because 
in this case, the leading-order electron-hole pair creation processes are active. 

Finally, we note that for even lower bias voltages, |e$| < 2ei, the current and the respec- 
tive current-induced vibrational excitation vanish in the limit A — >■ 0. 
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3. Resonant absorption processes mediated by a higher-lying electronic state 



In this section, we address the question whether resonant absorption processes with re- 
spect to a higher-lying electronic state provide a cooling mechanism (29] that prevents the 
vibrational instability we observed for a single electronic state (Sec. IIII B 1ft . Thereby, we 
assume the energy e 2 of the second electronic state to be larger than e% + Q. If the second 
electronic state would be located within the bias window, |e 2 | < |ei + £l\ < e$, we would 
obtain the level of vibrational excitation, which results from transport through this state 
only, because the other electronic state decouples from the vibrational mode in the limit 
A -> 0. 

Model calculations shown by the dashed blue line in Fig. [3] support this conjecture. The 
current-induced vibrational excitation of a vibrational mode that is coupled to a lower- and 
a higher-lying electronic state remains finite in the limit A — > 0. The parameters of these 
calculation are the same as those for the black line considered in Sec. IIII B 11 but include 
a higher lying-electronic state at e 2 = 0.8 eV that is coupled to the vibrational mode with 
a coupling strength of A2 = — 0.06eV and to the leads in the same way as the lower-lying 
state. 

This behavior can be rationalized by the master equation, Eq. (jSJ). To zeroth order in A, 
we obtain for bias voltages 2e 2 > e$ > 2fa + O) 

= FR,nPib"io - Pooroo r L,ii ~ Poofoo ^ ^,22^2,^ 1 2 fx fa + ^(^3 - ^1)), (31) 

Ku 3 

= r L ,nPoo!oo ~ PwJo^R,n ~ Pwjo ^,22^2,^3 1 2 fa + ^(^3 - fi)), 

Kus 

= ^2^K,22X 2 ^ lly3 X2 ;U3V1 Poo^ f K fa + ^fa - ^3)); 

= J]r i ,, 22 X 2 t i , ily3 X V 3, 1 p^ ^(e 2 + fi(z/ 1 -z/3)), 

Kv 3 

using the basis functions |00)|z/i), |10)|z/i), |01)|i/i) and respectively. Thereby, the 

populations P0101 an d Pun are treated as second order contributions < A 2 , because the 
population of the higher-lying state requires a preceding resonant emission process with 
respect to the lower-lying state (~ A 2 ). Since in the latter two of Eqs. (EJTj) all terms are 
either positive or zero, these equations can only be fulfilled, if the populations Poqoo an d 
P1010 van ish for values of u where fx fa + $l(y ~ ^3)) 7^ 0. Thus, vibrational levels with a 
quantum number larger than fa — p L / R )/fi are not populated in the limit A — > 0. Thus, 
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vibrational excitation is finite in the limit A — > 0, if a second higher-lying electronic state 
couples to the vibrational mode with a finite coupling strength A2. 



IV. CONCLUSION 

In this article we have studied various limits of resonant electron transport through single 
molecules. Thereby, we have focused on limits that, due to electronic- vibrational coupling 
of the molecular bridge, exhibit vibrational instabilities, i.e. an infinite level of vibrational 
excitation. To this end, we have employed a master equation approach that is based on a 
second order expansion in the molecule-lead coupling. To the given order in the molecule- 
lead coupling, this approach treats electronic-vibrational coupling exactly. 

An infinite level of vibrational excitation is trivially obtained in the high-temperature 
limit, k-gT — > 00, and/or in the static limit, Q — > 0, as energy-exchange processes transfer 
the thermal excitation of the leads to the vibrational degrees of freedom of the molecular 
bridge. A less trivial case is the limit of a large bias voltage, for which we have shown that 
vibrational excitation diverges with increasing bias voltage, $ — » 00. We have, furthermore, 
pointed out that this phenomenon results from the suppression of electron-hole pair creation 
processes. 

Similarly, electron-hole pair creation processes are blocked in the limit of vanishing 
electronic-vibrational coupling, A — > 0, if the bias voltage is large enough (e<3>/2 > e\ + fi), 
and consequently, vibrational excitation diverges in this limit. For lower bias voltage, where 
the leading electron-hole pair creation processes are not suppressed, we find a finite but 
non- vanishing level of vibrational excitation in the limit A — > 0. In the resonant transport 
regime, the vibrational degrees of a molecular junction are thus always excited, even if the 
electronic-vibrational coupling A becomes very weak. Cooling mechanisms induced by a 
higher-lying electronic state, as already pointed out in Ref. [291 ] . or coupling of the vibra- 
tional degrees of freedom to a thermal bath, as pointed out in Ref. 28j, however, prevent 
the vibrational instability in the limit A — > 0. 

Our analysis based on a master equation approach includes no higher order processes, 
such as co-tunneling processes. In the high-bias regime, such processes provide an equal 
number of additional excitation and deexcitation processes. Therefore, co-tunneling does 
not affect the corresponding vibrational instability. For a finite bias voltage, however, off- 
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resonant electron-hole pair creation processes provide an additional cooling mechanism. The 
role of these processes will be the subject of future research. 
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